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ABSTRACT 

A general survey is undertaken of the techniques used 
in determining loss functions to be minimized for the purpose 
of finding optimal (s,S) values in the Arrow, Harris, and 
Marschak dynanic inventory modes. A brief discussion of the 
cost functions and the aspects of demand, backlogging, lag 
time, and discount rate are presented in the interest of 
better understanding of their use in the real world, Analytic 
gOlutions to the model have been derived by two different, 
but related, methods involving stationary costs. The first 
method involves direct use of the cost functions in a 
Markov process to arrive at an integral equation of renewal 
which is developed into a stationary loss function by an 
Abelian limit theorem. The second method involves the 
determination of the stationary distribution of stock level 
through renewal theory ; the loss function for a represent- 
ative period is brought into the computations independently 
in determining the expected loss with respect to the 
stationary distribution of stock level. The two methods 
produce the same result. 

A digital computer simulation of the Arrow, Harris, 
and Marschak dynamic model has been outlined by Mo Geisler 
[ 2 ]- Some highlights of Geisler's work have been reviewed 
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PREFACE 


Inventory control organizations of the Navy have cone 
tracted for numerous studies into the decision problems 
involved in dynamic inventory for the multi-echelon systems 
a= found in the Navy. Also, a contract for the construction 
of a computerized dynanic model was first let in early 1955 
to the Stanford Research Institute. This work continues 
today with improvements being added asbreak-=throughs develop 
in the study of the nature of complex Navy inventory opera- 
tions. Industry has no parallel in complexity and magnitude, 
Models in use prior to the early part of 1955 were more 
applicable to reasonable-size firms. Consequently, the 
DbDoundaries of development of inventory theory are now being 
pushed by the military organization in keeping with the 
need for opener efficiency and economy. The magnitude of 
the difficulty in finding models that closely approximate 
the military inventory management problem are appreciated 
by relatively few. The optimization of buying, allocating, 
and pedi strieaedne based on complex statistical fluctuations 
in demand, order-arrival, and military loss due to shortage 

in a multi-echelon structure of centers and depots is an 
almost insurmountable obstacle to solution for a large 
part of the total carried inventory. 

The /exposition presented in this paper has been kept 


ji} 
Beistivtly Aumpie. Its purpose is to illuminate the basic 





concepts of the most simple dynamic inventory problem as 
rirst 0 Pare by Arrow, Harris, and Marschak and, thereby, 
gain sufficient insieht for proceeding into nore complex 
dynamic models. 

Source material for this paper has consisted primarily 
of references [2] | 2 lala and notes taken during a course 
of instruction in logistics given by Professor Thomas E,. 
Oberbeck at the United States Naval Postgraduates School. 

I wisn to express my gratitude to Professor Coerbeck for 
nis support and encouragement and to mv wife, Helen, for 


ner loyal devotion to the typing and clerical taskso 
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TABLE OF SYMBOLS 


(Listed in the order of their use in the text) 


reorder point of the stock level 

upper limit of the stock level 

stock level at the beginning of period t 
ordering cost per unit of item ordered 
reorder quantity in general 


loss or cost function for the nth 


period 

stationary or equilibrium cost function 

same as Ly) without regard to time 

stationary density function of stock level 

density function of demand 

stock level at beginning of period t after reording 
reorder quantity at beginning of period t 

penalty cost per unit of item not in stock 

holding cost per unit of item in storage 


cumulative distribution function of demand 


minimum discounted expected loss over an infinite 
time period with initial stock level, A ae 
discount rate which discounts future losses to 
present value 


nefold convolution for distribution of demand 


a random variable having nonenegative integer 
values which specify the number of periods less 
one in a reorder cycle. 


random variable of the excess distribution or 
the amount by which stock level goes below 
reorder point s prior to reordering 
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Onn ) @ random variable representing the sums of 
J demand just prior to D which causes 
a reordero v(y)+1 


Pr{ Fs-s>u| probability that excess of quantity of stock 


below the reorder point s is greater than 
some value Ue 
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foe labrodcuction 

The history and general characteristics of the inven- 
mery provlen under various operating conditions have been 
adequately discussed in pn. The (s, oo policy is widely 
used in practice and has been subject to extensive analyti- 
cal study. This ordering rule is simple and has the follow- 
ing advantages: (1) considerable intuitive appeal, (ii) large 
anounts of data available, (iii) a certain ease of proba- 
bilistic conputation and analysis and (iv) existence of 
analytic proofs as to its optinality under a wide range of 
conditions. Note that there are two approaches to the inven- 
tory problem. One is to choose a sinple policy (e.g. (s,S) 
type) and proceed with the analysis of specific (s,S) policies, 
amo the other @pproach is solve the entation problem through 
deternination of the characteristic of the optimal ordering 
policy which. minimizes or naximizes an appropriate objective 
function. Finding the characteristic of the optimal policy is 
very difficult. Thus the usual approach is to choose the 
mractical endl seu (s,S) policy and then proceed to find 
been Chis policy class a specific policy which nininizes 
the expected long run costs (i.e., specific values of s andS). 
Then the study of the inventory model becomes a study of the 
associated stochastic process which is generally a Markov 
process where the states of the system (current stock level) 


depend on the past only through the present. 


loften called the two-bin policy. 
1 


The Arrow, Harris, and Marshak dynamic inventory nodel 
(hereafter called the AHM Model) is the generic nane often 
cpiven in the literature to identify the (s,S) policy nodel 
where ordering takes place at the beginning of the period 
maen the stock level is less than the level s (i.e. y, <a) 
and demand is identically and independently distributed in 
each tine period. There nay or nay not be a specified 
fetivery lag. In this thesis, deliveries will be assumed 
to arrive innediately after ceing ordered. It is assuned 
that the reader is faniliar with the AHM dynamic model of 
mererence [ 5]. mone Simple restrictions on the cost 7uac= 
tions (having to do with the convexity of the expected values 
of the cost functions) according to Scarf [20] are the only 
Bewrictions required in order that the optinal poliveyeior 
the AHM dynanic nodel be of the (s,S) type. Thus over a 
wide rangs of conditions (including any demand distribution 
functions), there can be confidence in using (s,S) policies. 
Sections 2 and 3 are devoted to a resume of Cla lod enolic Funct aes 
Sroeother input factors in order to develop a vdetter under= 
Sbending of their neanine in analysis ~-and their relationship 


to the real world. (Pes 


| 





ee eesh = onct LONSs. 
The inventory problen has three cost functions which 
can have four forns (linear, convex, concave, or nixed), 
a Ordering Cost. 
puUppOSe 2 Represents the amoumteer stock to be 
ordered and c tne cost per unit ordered. A linear 
ereering= cost TSed rectly proport lomaeror je snout 


ordered, i.e. c*z. The case where c(z) is concave 
is a conmon situation in buying since this néans 
Ghat “each ecditiowal item costs less. Avspecia ae 
of concavity exists wnen c(z) is conposed of cez+tK 
where K is the administrative costs of processing an 
ordér. The convex ¢&se is unlikely to occur since 
it neans that each additional iten costs nore. The 
mixed case where pieces are nade up of linear and 
Sencayée porpTens is practical. 
me oOLOraeenvoat 

Usually storage costs are proportional to the 
size of stock on hand, but the cost nay increase nore 
rapidly if the storehouse capacity is exceeded and 
additional space must be rented. Storage costs may 
include such things as stock maintenance, storare 
ental, stock obsolescence orpyspolilage, stock repair, wee 
Gn Prenaity Cost. 


These .costs arise when demand exceeds supply 





as opposed to storage costs which arise due to 

supply exceeding demand. The forn of the penalty 

cost function varies widely since so many conditions 
exist. The sinplest is the linear case. A realistic 
case for industry is the convex function since a 

small shortage has snall consequences but a large 
shortage produces increasingly preater difficulties 
for the custoners. In sone nilitary applications 
Corie tae ibens Cirect is aifect the LS atoms short- 
ace of these itens will cause failure of the nission; 
consequently the shortase function could be considered 
a positive constant when a shortare develops and zero 
otherwise. In reality, the military supply items vary 
in importance and end-use. This is why a military 
supply system is hardput to define the nilitary essen- 
tiality of all items. So far, the military essential- 
ity designation has been applied to the allowance lists 
of designated submarines; however there is a nove afoot 
to expand this progran to include all naval vessels. 
The obvious difficulty in this project lies in review- 
ing and updating the allowance item essentiality 


designation. 





3. Other Faetors. \ 

Inventories are held for the ultimate purpose of 
satisfying demand, and money laid out for orders and items 
of inventory will not be returned for awhile. The alterna- 
tive to having money tied up in inventory is to have it 
invested in securities, bonds, or similar instruments, 
Long=term government bonds are a good standard for 
comparison since they are both secure and yield a reasone- 
able return of %. Lag time is the time between ordering 
and receiving, and it directly affects the amount of stock 
to be carried. Backlogging is the procedure whereby demands 
in excess of stock (unfilled orders) are kept on the books 
until they can be satisfied thereby giving meaning to 
negative inventory values, 

ae Discount Rate, 

A return of $1.0) at the end of a year is equiv- 
alent to $1.00 invested today in inventory. The discount 
rate is (Y= 1/1.0h which is an equivalent way of consider- 
ing the comparative return. Thus a $1.00 return in a year 
is equivalent to ( dollars invested in inventory today. 
Considering that the $1.00 return on investment is re- 
invested each year or remains invested at %, then XX, 
as, ax, err ewvers pepresents the discount rate for 
$1.00 after one year, two years, three years, etCo he 
(Loe. a® is the present value of #1.00 which has | 


{ 
been tied=up in inventory during the nth period). =: a 





ror all practical purposes (Y= 1 when the time horizon 
is short. 
be Demand, 

Demand is usually assumed to be independent 
of the decision makers control; however, this is not 
always the case in industry (e.g. demand is directly 
related to advertising), In military inventory it seems 
safe to assume that demand is independent of the decision 
maker. Demand may be regarded as deterministic in future 
time periods. However, it is more realistic to regard 
demand as probabilistic. Although the demand probability 
distribution may be known to change such as in seasonal 
fluctuations or long-term trends, the state of methodology 
imestochastic processes for inventory control requires 
the assumption that demand is identically distributed 
mmc me time periods. EFecause of this restriction, the 
deterministic case is sonetimes a better approximation 
for demand when it is rapidly changing over time (i.e. the 
deterministic case fixes demand in the various time 
Periods according to some previous knowledge of its 
behavior) 6 

Another way of looking at the aforementioned 
probabilistic case is to allow the size of each demand 
to be a random variables Consequently, this is a = 
tinuous time Stochastic process which is more appropriate 


when orders can be placed at any moment of time. For 





application of this concept, the assumption is made that 
nea are independent and identically distributed random 
variables. 

When nothing is known of demand one may apply 
decision-making under uncertainty as developed by Wald 
[3] which involves both estimation and decision-making 
as successive observations of demand are taken, a com- 
plicated sequential decision problem. Y. Fukuda provides 
a qualitative approach to this solution [4]. 

Co Lag Time. 

There are three cases of lag time to consider. 
fim this these the time between order and delivery is 
assumed to be zero in order not to over-complicate the 
computational techniques to be displayed; however the 
extension of this model to the case of fixed lag time 
between order and delivery is presented by Karlin and 
Scarf [1] and has the same characteristics as the dynamic 
inventory model with no lag time when both models assume 
backlogging (see subsection 3d). The case where lag time 
is a random variable ee been treated by Scarf [1] (see 
emeerer 16). 

It should be noted that when demand is deter- 
ministic, the consideration a fixed lag time is in- 
consequential since we Sviny order sufficiently in 


advance to cover the known demands. In all cases the 





lag time directly affects the amount of stock held 
(i.e. an increase in inventory to protect against shorte 
ages until deliveries arrive). There is onepractical 
consideration in the dynamic model where lag time is 
a random variable for which no solution has been found 
through the dynamic programming formulation at the 
present time and that is the condition where orders 
arrive out of sequence. This is referred-to as cross 
over. Remember that order quantity is a random variable. 
do Backlogging». 

When demand for an item exceeds supply, then 
the policymaker may do one of several things: (i) make 
a premiun order thereby taking a penalty in the amount 
of the premium costs, (i141) do nothing and allow the 
customer to go elsewhere for satisfaction thereby taking 
a penalty in customer goodwill, (iii) backlog the order 
awaiting arrival of orders to satisfy the unfilled demand 
thereby taking a penalty that may relate to both customer 
goodwill and/or increased ordering costs. In this case 


a negative inventory becomes meaningful. 


ho Methods of Solution Available. 

At this juncture there exists a stochastic process 
based on the (s,S) type policy which may now be treated 
in one of two waySo 

ae First Way. 

The first way consists of viewing the transient 
behavior of the stochastic process through neperiods 
or an infinite number of periods. This essentially is the 
original approach used by Arrow, Harris, and Marschak 
[5], [>] in setting up the loss function for a general 
period t, then sunming the expected losses over a given 
time horizon with respect to the random variable of 
stock level, The distPibution of stcek level is deduced 
through a functional relation with the random variable 
of demand. This method will be outlined in section 5 where 
it will be pushed to an approximate solution through use 
of an integral equation of renewal theory [2], [5]. 

De Second Waye 

The second method of treatment deals with the 
stationary phenomena of the process. If yy, represents 
the stock level at time n, then subject to mild stabile 
ity requirements it can be shown [1 (refer to.. pages 
234-237 and pages 292-297) that y, possesses the so- 
called ergodic property (i.e. the distribution function 
for y, converges to a limiting distribution function) 
when the optimal policy is of the (s,S) type. Thus the 
limiting or stationary distribution is a function of 


o 





the policy and the demand distribution, and it is in- 
dependent of the initial stock level. The beauty of this 
method is that the cost functions have not entered into 
the problem up to thik point. Using the loss function 
representation for a general period n, denoted Liy.)s 
and dropping the dependence on n, the expected value 

of Py) with respect to the stationary density yields 


the equilibrium average loss, denoted as follows: 


Lam | Lo g(y) dy 


where gly) denotes the stationary distribution for the 
stock level. | 

The stationary distributions discussed in subsection 
eb can be found through the application of renewal theory 
[2] (refer to, chapter 15) or through the application of 
differential equations [2 |(reter’ to, chapter 1) when 
the demand distribution function is a member of the 


Gamma family of distributions, 


: “TR 
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ik Original Method of Solutione 

It is proposed that the original Arrow, Harris, 
and Marschak dynamic model be solved for the optimal 
values of (s,S) where the demand density function is 
specified as £(€) =@E™> . The detailed computations 
are straightforward and can be found in Appendix I. The 
rule of action for time period t with Y, a8 the steck 
level, xX, as the stock level including replenishment, 
and ey, as the demand are represented both analytically 


and graphically as follows: 


If y,>s, then 2,=0 (and X,=Y,) 


Ig Ye=s> then 2,=SeJy (and x,=S) 


X5 Xx 


td 
Wr 





time 
Steck Level\as a Function of Time 


Figure 1 


Let c(°), p(*), and h(°e) Ls the ordering, penalty, and 


holding cost functions respectively. This model has cest 
| 
pel 
a 
} 





functions defined as follows: 


c(e) = K 

. 
p(e) =p 
h{e) = hxy, 


Under an (s,S) policy with the initial stock position, 
Je? specified, there is a Markov process on the sequence 
of values for the stock level. The loss, Leyy)s in any 
peried (t, t#1) depends on the stock level, Ths at the 
beginning of the period, the cost functions, and the : 
demand distribution function, F( € ye 
(501) Liye 7 * B-Pzg1). tom amas 
hxS + p [ 1-F(s)| for y,<8 
An important feature of the function ey), which will- 
be used later in finding the optimal (s,S) values, is 
that the function involves s and S as parameters ana 
is constant with respect to Y_=53 thus Pto) - fis) = Ke 
Consider L(y.) as the minimum discounted expected loss 


which will be incurred during an infinite time period if 


} 
ce is the initial stock level and an optimal ordering rule | 


: 
is used throughout this infinite time period. If an order » 
at the first stage brings the stock level up to amount 

; i 
y and an optimal ordering policy is followed in /the 


second stage onward, then the expected loss from’ the 


second stage onward discounted to the present iss 


12 


hy 


S 
‘oe [ L(s- €) aP( £ )41(0)(2-"(5)] if y<s 


O- 


eaec) 


bf 
X [ L(y- &) aF( € )+1(0)] 1-P(y)] if y>s 
0- 
Censequently, the minimum expected loss is achieved by 
the policy that minimizes the sum of the expected loss 
for the first period and the discounted minimum expected 


losses in the future: 


S 
ee) Ll y)= Linea | L(s-&) ar( £)+OL(0)[1-F(s)] 
O- 


Lis 
| sy 
(Soh) Lins Lenear| Ly-&) aP( € )+Q1(0) [1-F(y)] 
0 


if y>s 


These are the same functional equations (hell) and (hel.) 
of [5 | showing the recursive properties of the Sygate 
model. The lower limit of integration is set at O- in 
order to insure that the random variable for demand, é, 


4 


5 
‘does not have a positive probability at 6 =o, and, there- 


by, avoid a discontinuity in the Stieltjes integral . 


Ny 


' af 
The exposition in this section could be shown for the \ 


case in which the demand distribution is discrete and 


all functions are defined on integer values; the intearetey 


would be replaced by appropriate summations. . 


13 ( 


J 


t 
\ 


} 


/ 





It was noted earlier in equation (5.1) that Ly) 
is independent of y for O#y£s. Note in equation 
(5.3) that L(y) is independent of y fer 0Hy=—s,. 
Putting y=0 in (5.3) and y=S in (5.4), then subtracting 


the latter from the former produces: 
(5.5) Lio) - L(s) = fo) - £(s) = xk 


Equation (5.4) can be solved for the function L(y) by 
eonsidering L(0) as an unknown parameter, then (5.5) 
can be used to find L(0). But before this solving takes 
place, (5.1) is changed into the usual form of the 
integral equation of renewal theory whose solution is 
found in terms of L and f. It is ae solution that is 
used to find L(0). Substitutions, change of variable, and 
renewal summation of convolutions for demand provide the 
means for finding a solution [ 5]. 

The steps begin — breaking-yup the integral on the 
right side of (5Solt): 


y ys y 
[ L(y- &) sf )=f L(y-€) aF(€ +f L(y-& )aF(é ) 
O 0 


yes 
ys i 
-| ur Ge € rsuco) | dF ( & ) 
0 y-s 


In the last term L(0O) = L(y-s). ‘Substitute this 


Ly 








inte (5.4), then in the next step substitute 7)= y-s and 
L(y) = L( 7] +s) =)(N)): 


| — y 
io.) ugieLisr+0xf Hiy- §ert@eaxnion] ar(€)+CxL(0)| i-F(y)] 
0 a 


" 
AMP =LN+s) + A(N-§) aF (E)+QL(0) rty)-r0pe[a-ets)} 
O 


n | 
AMpehnes+a} AN-Ear(E)+axL(0) [1-F(N)] cer D0 
0 


This is the integral equation of renewal theory. The n-fold 
convolution of the distribution of demand is F (€) which 


is used in the solution of (5.4) as follows: 


CO. 
(5.6) HOE) =) om AE) 
n=l 


where the convolutions are defined as follows: 


¢ 
F(é) -[ F (eu) ar(u) for n=2, 3, * °°. 
0 


The solution te (5.4) in terms.of L and pis shown on 


page 33 of [| 5 Je 


LS 





= 


8 
K + As) -[ K(s- € ati ( €) 


0 
(1-0) [14H (s-s)] 





As pointed-out on page 3h of [Ss]. the determination of the 
optimal policy which minimizes L(0O) is tantamount to find= 
ing the optimal policy that minimizes L(Y¥o) provided that 
the optimal value of s is greater than zero. The finding 
of the optimal (s,S) _. required that the system of 
equations 2-110) = O and Sate = 0 be solved in terms 


- rae ; 
of S and A= Ses. Also <= O) must be positive for the 


minimum to exist. The appropriate equations are: 


ale 
” 


: S-s 
(5.9) Lo)=nepee(s)+] [nepxt (s-E aH (€) = 0 | 
0 


Q/|OV 
> 


S-s 
L(0)=h(S+s)=p| f(S)-F(s)| vce] [ h=pxf (s-€)} Hoag 
0 


Tt is assumed that H a §) fs a '*? BiG. ie 


One can find hy () for many demand density functions 
which have practical application to actual demand behav- 
for, However, solving (5.9) for the optimal (s,S) values 
may require numerical methods. 


The other alternative to this method is a simulation 


: 16 





ee) 
of (5.8) using predetermined functions of dH ay (€) =) XR (E)« 
For some guidelines on the number of runs 4 Oe ae a 
95% confidence level, see Geisler [2 |. An attempt was made 
to setvup this simulation on the CDC 160) computer using 
FORTRAN machine language; however, unfamiliarity with 
FORTRAN caused the effort to be time consuming both in 
preparation of the program and in debugging. It is suggested 
that the next undertaking of this problem begin very early 
in the learning phase of programming and that another 
machine language, such as SCRAP, be investigated for better 
adaptability of this problem. 

As foretold in section 1 the equation is now trans- 
formed to a stationary cost: equation which implies ene 
meemitoss in period t, f.. approaches a limiting value e. 
pe. as too and is independent of e's the initial stock 
Orel. By use of a standard Abelian theorem | 20 | and as 
heuristically outlined on pages 37-38 of [S| the right 
side of equation (5.8) may be multiplied by (1-C%), and 
letting X~»1 as suggested by fin (l=) L(0) =P oo the 


stationary loss equation is obtained: 
S-s 


K + Xs) | Pis-&) an(€) 
0 


1 + H(Ses) 


! 


Mee ~ 





A solution using demand density function ral ) = rs wie 


Ii 
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1 ual 





be found so that a comparison can be made with the 
solution in section 6. On page 39 of ie. demand dene 


sities of the type 


k 


« § 
f(g) A fabie? du, k>0O and B>o 
0 


(where k is an integer) are shown to give the following 


results; 





1/k {/k 
eae al"? Lida 7 ls & Bu 7 

OL . j 

O 

The gamma type demsity function is chosen since the product 
of the choheet eae function netimes yields another 
gamna density with parameters nk and B (the sum or 
addition of n independent identical random variables 
is the product of the characteristic function n-times.). Here 
WI, , e/,, seeeeereo y (po) are the yeh 
For Bie ) = aa » the parameters are k=l and [Se 


(thus (g/ 


roots of unity. 


1 1). Since M=1, then 


and thus, 








As’ shown in Appendix I for fixed A= S-s, the minimum 


occurs at 


(! 
li 


fa, — = Meg) + A 


s = dn — - Daiis A) 


However, it is useful for comparison with other methods 
of solution to find the optimal value of A. From 


Appendix i 


2K 
A =Yy— 
4 
_g hh +Y2Kh 
Ee = 
p 


which will turn cut to be the same solution in section 6 
wherein tne stationary distribution for stock level was 


UsSede 


AES. 





6. Selution Using Stationary Distributions of Stock Level. 
As noted in chapter 15 of [ 1], Looe Git (1-00) L(y) 
may be also determined by finding the stationary distri- 
bution of the stock level and afterwards calculating the 
costs for a single stage problem where the stock level 
is a random variable having the stationary distribution. 
The stationary distribution may be obtained for a general 
demand density function f ( €) and is composed of two 


sections. For stock level, y, in the range [s,s], 
iors 1.) g(y) =C 7 Baio y ) for s#yeés 


The constant C is determined by the appropriate boundary 
conditions and BO is the n-fold convolution of the demand | 
distribution function (cf. equation 5.6). The stationary 
density function for y<s is a multiple of the so-called 
excess distribution for the renewal process of independent 
random variables ae Although (6.1) may be difficult to 
find, the distribution for the excess variable is more 
readily obtained by use of Laplace transforms. The 

excess variable can be characterized as 
Se + ¢, eicieiais|s co oeiien Sein ~ (Seg) oa AEP 
where v(y) is a random variable whose falue may be an 


integer from the set (0,1,2,3, eee+++) satisfying the : 
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follewing inequalities: 


\ f +g, + eoecces oi) 


¢ +¢. + eeeeeo +o yy41 = S8 


Letting 


: P 


“Ss = a ce + ceeeee re and 


y) 


Pryiyyar 7 $1 : ae Sv(y) 41 


the excess variable of y may be pictured as follows: 


1 
ere 
Ses : 
Pry) Dycy)a1 
Excess Variable of Ses 


Figure 2 


The renewal type equation follows from this definition of 


ao 


Pr 7.>+]* | Pr To gay 2 eh (a) [2-F(s-s-u)} 


A renewal process is defined as a stochastic process 
compesed of sums of positive identically distributed and 


independent random variables, This much of the theory 


a 





is presented in order to provide sone intuitive appeal. 
Those interested in the detailed study of renewal processes 
are referred to chapter 15 of [a]. 

When the demand density functions are members of the 
gamma family, the stationary density function of stock level 
may be derived alternately through the use of differential 
equations as in chapter 1h of [2]. The stationary distri- 
bution, gly), for the case where f(é)= @€ “is taken from [2 | 


and is found to be 





for s<y=S 


er (s-y) 


aie for ys 
lS G=ns 


The stationary loss a a straight forward calculation 


Le -f hove (y) ay 


where the one period loss in our case is set-down in (5.1) as 


Li = hy + p[1-F(y)] for s<y=S 
—_ hS + p|1-F(s)] for yss 
Since F(€) = 1-@ 5 then 
in hy + [er | for sey<s 
hS + p[ e + K for ys 


Sinde ‘our model allows backlogging, negative stock levels 


ee 





may exist, and then 





S S 
7 1 eae a 1 
ian | (hsp eae e 3 Page | (ny +P e 7) 
- S 


a 1+S-s +Ses 


K+hS + pé + 3n(S-8*) 


1+Se-s 


This is the same stationary cost function as (5.10), and 


the solutions are the samé. 
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fer The Simulation Problem. 

The great advantage in the use of renewal theory is 
that the cost functions enter separately into the solution 
and in addition they enter into a simple expected value 
form with respect to the stationary density of stock level. 
Thus the costi functions nay be more complicated without 
unduly icden tae the solution. This approach also provides 
a more flexible nethod for simulation wherein the program 
will not be ‘tied down to a specific set of cost functions 
as would te the case in a simulation of the equation for 
lO), Painted cut earlier in section 5. 

Geislerts ie cort [2] provides a statistical evaluation 

| 


of the sample Te required for a desired level of precision 


Lor estimating: hortages and overages (i.e. stock level 


above zero) in'the AHM dynamic model for both lag and ne 

. lag in delivery with linear penalty and holding cost and or- 
eee cost, as em The precision represents a way of 
finding sample size when the eee of shortages or 
overages is within K% of the true value with 95% confidence. 
Part D of Geisler's report has sections iitiecend iV watch 
analyze the given nodel by invMeEeute once and by time 
periods,respectively. The tide pap iedumode1 is what has 
been discussed in ‘ie este: The inventory cycle is the 


number of periods between reorder where the number of 
? : 


periods is a random variable. The concepts in renewal theory 


2h 





are covered in chapter 15 of Fue sections III and IV. 
Sections III and IV outline the arelytic solutions. 
Subsection III-3.b shows the specific expected costs with 
respect to the stetionary density of stock level when 
oo i) =\e e« In subsection IV-2, the expected values of 
shortage level, ston level (overage), and reorder quanti- 
ty per period are given on pages 79, 85, and 89 respec- 
tively. To get the stationary loss (cost) nultiply, 
meeroectively by po, h, and c. To the reorder cost 1c 


must be added Consequently Geisler's report has 


toa 
xO FHA) 
me 7.ded an exercise reference in the study of renewal 
theory and its utilization in finding an analytic solution 
to the AHM dynamic model with less restrictive cost functions. 
aeemmeone stationary costs have ceen derived, either directiy 
Sumeyeuse of numerical nethods, a Ssinulation can be under-= 
Geren to find the optinal (s,S) policy which cinta the 
Stationary loss function. Section V of part D can be used 

\ 


in the simulation to deternine sample sizes. 
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APPENDIX I 


9. Computations for the Original AHM Solution of Section 5. 
Starting with the case where (f)=er$ , the values 
of (s,S) in functional form will be found that minimize 


K+hSe p@n~® +2n(s°-s2) 


CoS 





1l+Ses 
Let / =Ses. 
(1+S-s)(h+hS)- [Kenstp @7%e 3n(s°-s?)] 





oe 
® 35 Leo (19S=s) 

+ (145-8) (-p@*-hs)=[K+nS+p@ +4n(S2-s%)] (-1) 
© 2 f= 2 

os (1+S-s) 


@D' hthAtns+ns A-K-hS-per® -An(S+s) = 0 
hth /A\+3Ah(S-s)-K-per =0 


~K-p@ S+h+A h+3 Ah = 0 
i 


@' -pe@*-Ape *-ns-A nsdK+hStp@ *+3Ans+d Ans = ( 
-Ape +h(s-s)+eAns-3Ahs+K = 0 
-K = -Ap@ “thA+3Ah(Ses) 
-K = -Ap@ +hA+inA® 
“Substituting ©' into @' for -K: 
(A p@ *+hA +bh A“) -p@S+hth A +3H ARM =O | 
~p@ "(1+ A)+h( A“+2A +1) =0 | 
pe (1+A) = hla+A)? 


at 








MD 
u 


h(1+/A) 
st = fn — = Lr (1+ A) 


s* and S* denote 


si all : Dn nee iu optimal values 

h 
These results agree with page hO of [5] where A has a 
fixed value. In order to solve explicitly for A , @ 
mathematical device must be used. This device consists 
of removing the variable h by taking the less function 
with respect to h after certain substitutions into Wea 
Substitute the following into Loo? 


pe = h(1+ A*) 








; p 
3s* = 
h(1+ A*) 
- Pp ‘ 
S* = ANG 
h(1+ A*) 
p : p 
K+h| fv ——_——— + *]+n(a+ At) +b A* 2 la —_——3 7 
| h(1+ A*) A h(1+ A*) 


= 





1+ A* 


“ p 
K+h E A*t+2 Ara +h(1+ We) fe _—— 
f oui. 4 
CO = 





i+ 7 aNes 
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2 2 
K" x \it2 el 
foo _ yn = + 1lix———_-—- + fn p" = fa(1+ A*) 
h i+ A\* 1+A 





where K" =— and p" =— . A check of the consistency in the 
h 


units of the equation shows that 1! = (one charge-unit/period) 
= (storage cost/unit-period) (no. of units stored, averaged) o 
In the first term the order handling charge is divided 

equally among the number of units in the dynamic range of 

+1. The third term is the penalty premium expressed in 

(1' units x1 Unit) and averaged. This analysis of consist- 
ency helps verify the propriety of the substitutions. Now 


2 
minimize with respect to A : 


3 g -K" (1+ Ko (Zee ee At +2 A*+#1) Lie 
Ba ee ee ee 
oh as A? ana A? 


é 
~K"4 [A *¥ +3 A* +2 -2 A**-2A*-1]-1 +A* = 0 
-K" + BA = 0 


ek 
A” = VaR -y— 
h 


To set the (as) values after A* is determined, an 


equations for s* is used: 


pe” = h(1+A*) 


Lt fi OK h+ ¥2Kh 
e =— (1+4— ) = 
p n P 


This completes the analysis. 
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